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1 Weak Values and Variances 

1.1 Weak Values 

1.1.1 Introduction 

In the two state vector formalism, a weak value is equal to the transition 
amplitude between a prepared initial state and a postselection state [3]. We 
focus in this study on the weak value and the “weak variance” of the mo¬ 
mentum operator p postselected by position (ccj, 

{x\pH) 

xPl/i - / I /\ ’ xPljl G 

A number of the following results should generalise to any pair of canonically 
noncommuting variables, although a general definition of the weak variance 
is not the objective of this study. The weak value xP'ip G C can be split into 
real and imaginary parts [4] as 

One might naively attempt to define a weak variance as the second central 
moment out of weak values xP^ — xP^j'^- However, the fact that weak values 
are complex requires one to take real parts of this naive construction - and it 
makes a difference whether one squares the second term before or after taking 
real parts: Re{xPjp^) ^ {^^xPipY- This ambiguity motivates the search for a 
more ‘fundamental’ approach to these statistical quantities, which this first 
section will attempt to describe. 

This study is organised as follows. We begin by reviewing the relation 
between the real parts of weak values and the phase space formalism of 
quantum mechanics. The quasistatistical interpretation of this formalism sug¬ 
gests a natural and unique definition for the “weak variance” in terms of the 
Wigner function. In §2, the embedding of configuration space formalisms into 
the phase space formalism then allows the weak variance to be related to the 
thermodynamics of the Madelung fluid, to the de Broglie-Bohm quantum po¬ 
tential Q, and to an experimentally measurable combination Re(a;P^ — xPip^) 
of weak values. 

In §3 we discuss the titular problem, that the negativity of the Wigner 
function can result in a negative weak variance. The weak variance there¬ 
fore does not define a standard deviation of weak measurements from the 
weak value. The sign of the weak variance is shown to be related to the local 
extrema of the probability amplitude. The relation between weak variances 
and determinism paves the way for a detailed study of the Classical Limit 
in §4. We derive the classical limit from a variational principle involving the 
average stochastic contribution to the total momentum variance, and demon¬ 
strate the existence of a semiclassical limit where particle motion occurs, on 
average, classically. 
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1.1.2 The importance of the measurement apparatus 


It is important to note that weak values are not measured independently 
of an interaction Hamiltonian or a meter system. The real part of a weak 
value is well known to be the conditional average one would measure in the 
ideal limit of zero disturbance form the measuring apparatus pimiH] . The 
imaginary part of a weak measurement is unrelated to observables, it arises 
from the disturbance due to coupling with a measurement apparatus [B]. As 
such, it is not uniquely defined and not intrinsic to the measured system [^. 

Since only the real part informs us about our operator as an observable, 
we opt to describe the variance of only the real part of the weak value in 
this study. However, the imaginary part in m will emerge naturally during 
our investigation of weak variances for the real part, without reference to 
a specific interaction Hamiltonian or meter system. In principle the ‘ideal’ 
decomposition © would be measured as [5115] 


xPip — 


V / 


+ r] Im 


f {x\p\f’) \ 

\ W) r 


( 2 ) 


such that the theoretically natural imaginary part in is arbitrarily scaled 
in experiments by adjusting the value of rj. This parameter depends on meter 
states and observables, as well as the interaction strength with the apparatus 

0 - 

We shall therefore discuss the imaginary part as if it were intrinsic to 
the system, absorbing all dependence on the measurement protocol into the 
parameter 77 . This has the disadvantage of illusorily disconnecting these re¬ 
sults from the experimental setups that give rise to them, which comes hand 
in hand with the hidden advantage of being a more universal treatment than 
that derived from any specific experimental setup. As such, the weak variance 
we define is related to the stochastic nature of quantum theory; it is not the 
experimental variance of any apparatus (although it may manifest itself as 
such in a measurement setting). Not without irony, this apparatus-agnostic 
treatment will suggest multiple possible experimental setups to study the 
weak variance f 112.31) . 


1.1.3 Weak Values as Conditional Expectations 

In the quasistatistical interpretation mm of the phase space formalism 
the canonical observables x,p gM. are treated as random vari¬ 
ables. Since these observable are incompatible, no joint probability distribu¬ 
tion can be reasonably be given for them - the phase space distributions are 
understood to be bivariate pseudoprobability distributions. For example, the 
Wigner function }V(x,p) is not positive definite. For a review of extended 
probability we refer the reader to m- 

Our choice of the Wigner function to define the weak variance is moti¬ 
vated by this quasidistribution’s special status amongst physically equivalent 
representations of the phase space formalism (which correspond to different 
choices of operator ordering, [x,p] 0). Indeed, only the Wigner function 

produces the correct marginals and observable averages by direct integration 
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[14) : other quasidistributions require explicit convolutions over noncommu¬ 
tativity terms. As such, the Wigner function is uniquely singled out by the 
quasistatistical outlook, where such noncommutativity terms would interfere 
with the interpretation of observable averages as moments of the pseudodis¬ 
tribution. 

The real part of a weak value is a conditional average IZIIHIIIS]. Hence, to 
obtain weak values of momentum, one must consider not the momentum p, 
but instead the mean value of p at a given position, for a statistical ensemble. 
One may consequently express the experimental protocol of postselection 
with conditional random variables: the quantity p\x {p conditioned by a:) is 
also a random variable. To find the distribution of p\x, we need to consider 
the conditional distribution 'W{jp\x). The Wigner function has a marginal 
distribution 

Wix) = J Wix,p)dp = \'ip{x)f , (3) 

and so the conditional distribution for the random variable p|a; is simply 


W {p\x) 


W (x,p) 

W{x) 


(4) 


These observables p, x are incompatible: in what sense can we condition 
one on the other if they are not meaningful simultaneously? A random vari¬ 
able is defined by its distribution. Individual realisations of x and p are mean¬ 
ingful positions and momenta only because their distributions (the marginals 
of the Wigner quasidistribution) are proper distributions. The random vari¬ 
able p\x is defined by the conditional quasidistribution W (p|a:), which from 
o is clearly not nonnegative. Hence individual realisations olp\x may occur 
with negative probability. This follows the same logic as stating that joint 
realisations {x,p) may occur with negative probability as realisations of the 
full Wigner function. 

However, the fact that individual realisations of p|x can occur with nega¬ 
tive probabilities does not preclude the use of this conditional variable as an 
intermediate step in a calculation of a physical quantity m- Furthermore, 
the mean or the variance of p|x need not be confined to intermediate steps in 
this way, since they are not random variables but instead contribute to the 
description of the physical state of an ensemble (i.e. they descibe the Wigner 
function as a whole rather than the outcome of a measurement). 

The mean value of p|x, 

Ep (p|x) = p, (5) 

can be thought of as a function of x cni- This mean value p will in ^1.2.11 be 
shown to equal the real part of the weak value xPtP ■ One may then express the 
conditional mean momentum in terms of W as the first conditional moment 
[T], i.e. the first moment of the conditional distribution: 


p = y pW (p|x) dp 


W{x) 


( 6 ) 


pW (x,p) dp. 


(7) 
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It is possible to formulate complex weak values of arbitrary observables 
directly in the phase space formalism using cross-Wigner functions m, but 
as noted earlier our focus is currently limited to the variance of the real part 
of position-postselected weak momentum measurements. 


1.2 Weak variance as the Second Conditional Cumulant 

Through the statistical Wigner-Moyal formalism described above and in [I] , 
the weak value p takes on the role of a statistical average. However, an 
average is scientifically worthless unless one can provide an associated spread. 
It is well established that the variance, i.e. the second central moment of 
a distribution, is equal to the second cumulant of that distribution. With 
some effort, the same can be shown of conditional variances and conditional 
cumulants [TH] , and extended to those of pseudodistributions like the Wigner 
function [T]. 


1.2.1 Derivation 


Using the characteristic function associated to W {x,p) and the polar decom¬ 
position tjj = the conditional cumulants of this distribution can be 

written for even cumulants (n = 2,4,...) as [T] 


'‘"'“^=(4) (s) 

and for odd cumulants (n = 1, 3,...) as 


( 8 ) 

(9) 


The derivation of this result has been reproduced in full as an Appendix. 

By definition, the first two conditional cumulants n = 1,2 are also the 
first two conditional moments, and one finds a simple relation [1], 

. d(hS) 

for the (n = 1) average p = Ep (p|x). Simple comparison of this expression 
with m should confirm that p = Re {xPijj)- The n = 2 expression then defines 
the weak variance of the random variable p|x: 

Vp(p|x) =Inp(x). (10) 


Dimensional analysis confirms that this quantity has units of [momentum]^, 
as expected of a (conditional) momentum variance. Because of the clean 
splitting between phase-dependent odd cumulants and amplitude-dependent 
even cumulants, the weak variance does not depend on the wavefunction 
phase, except indirectly (via the time evolution of the amplitude). 
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1.2.2 Other expressions for the Weak Variance 


This second conditional cumulant is also expressible directly as the second 
central conditional moment of the conditional distribution W (p|a;): 


Vp {p\x) 


J W ip\x) {p\x 


Ep {p\x)f dp. 


( 11 ) 


Furthermore, as we shall demonstrate in 112.1.21 the weak variance also 
takes the following form, more readily related to weak values: 


Vp {p\x) = -Re 


( {x\p‘^ |V>) 


V (^IV') )) ■ 


This last expression shows that despite our use of unphysical conditional 
random variables p\x in its derivation, the weak variance is indeed related 
to the intuitive form {p^) — (p)^ of a variance. It also shows that the weak 
variance is measurable in an experimental setting. 


1.3 The Uncertainty Principle 

It is possible to relate the conditional variance Vp (j)\x) to the unconditional 
variance V (p) that appears in the uncertainty relations. 

1.3.1 The Law of Total Variance 

One may readily verify that the law of total expectation, 

E (p) = Ea; (Ep (pIx)) = Ej, (p), (12) 

reads in our quantum setting 

{%j)\p\ij;) = jpdx. (13) 

In the same way that one derives the law of total expectation, one may 
derive the Law of Total Variance: 

V (p) = E, (Vp (p|a:)) + V, (Ep (p\x)) (14) 

= Ej, (Vp {p\x)) + {p) , 

where V^, (p) is the variance of the weak value p as a function of x, and where 
Ea; (V(pk)) is the average value of the weak variance (as a function of x). 
Yx (p) can be simply given by the marginal variance 

'^x(p)= J W{x){p-Ex{p)f dx, (15) 

where the law of iterated expectation o gives Ea; (p) = {tp\p\tf). Hence, 
Yx (p) measures the standard deviation of the weak value p from the usual 
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expectation value (p) as the postselection position x changes. Particularly, it 
does not measure the variance of p\x (or of any other random variable): the 
variable p in the variance Yx is deterministically defined by (1151) . The non- 
weak variance Yx {p) contributes to the total variance of p of the uncertainty 
principle, but it is not itself a manifestation of quantum randomness. 

By contrast the conditional variance, 

Yp {p\x) = y W (pIx) {p\x - Ep {p\x)f dp, 

is truly the variance of a random variable: it says how far from the mean 
Ep {p\x) = p one may expect to find the random variable p\x. Unlike the 
total variance or the variance of p, this “scedastic function” depends on the 
postselection x. It is in this sense that we identify it as the weak variance. 

These two parts of the total variance allow interesting limits where their 
contributions are negligible one with respect to the other. For example, in 
stationary states there is no variability in p{x) by construction, and the total 
variance of the state is supported entirely by the weak variance. A more 
detailed discussion of these limits, and their relevance to the classical limit, 
is postponed until the relevant properties of the weak variance Yp (p\x) are 
introduced ( 114.11) . 

1.3.2 Apparatus vs Fundamental Weak Variances 

This connection between the weak variance and the momentum variance re¬ 
confirms that the weak variance is central to understand the predictions of 
QM. Oddly, we find that the weak variance is related to the fundamental mo¬ 
mentum variance Y{p) of a state (as described in the Robertson uncertainty 
relations) rather than the momentum change Ap when a state is disturbed 
(as originally described by Heisenberg). 

Since Vp (p|x) is given by differentiation of the wavefunction amplitude 
(nni), it depends only on the values of the wavefunction amplitude in the 
neighbourhood of the postselection point x. This choice of postselection point 
is not apparatus-dependent, although issues such as pointer variance for the 
postselection variable will certainly affect the analysis of any experimental 
study. 

The weak variance is further distinguished from an ‘apparatus quantity’ 
by the fact that it is not an ‘observable’, in the sense that it does not corre¬ 
spond to a Hermitean operator (cf 112.1.d|) or equivalently (in the presenta¬ 
tion adopted above) to a real random variable. Instead, as a cumulant of the 
Wigner function, it describes the state of the system, irrespective of how it 
is measured. 




2 Connection to other formalisms 


One can recognise [2] the probability flux density f pW {x,p) dp = mj in 
and re-express this equation using the marginal , as 


P_ ^ j {x, t) 

m |^|2 


(16) 


The trajectories of the Hydrodynamical and Bohmian pictures are generated 
from (HU, by setting a velocity x = p/m and solving the resulting ODE for 
x{t). We will present only those elements of these formalisms that are of 
interest to our weak variance study; a more detailed account of this ‘embed¬ 
ding’ of configuration space descriptions of QM into phase space is given in 

m- 

The derivation of the real part of the momentum 


p = da; (hS) 


(17) 


from the Schrodinger Equation is present in almost every work related to the 
Hydrodynamical or de Broglie-Bohm formalisms since 1927 m, therefore no 
detail will be given here: The wavefunction is polar-decomposed, derivatives 
of the components split the Schrodinger equation into real and imaginary 
parts, from which (II3 is obtained by either appealing to an analogy between 
the resulting equations and the Euler Equations (in the hydrodynamical set¬ 
ting) or to an analogy between these equations and the Hamilton-Jacobi 
Equation (in the Bohmian setting). 


2.1 Bohmian approach to the Weak Variance 

2.1.1 Prelude 


In the quantum Hamilton-Jacobi formalism [201ICTI221I28] . Hamilton’s Princi¬ 
pal Function is S {x, t) = hS (x, f) -l-const., such that the polar decomposition 
(HZD becomes the Hamilton-Jacobi definition of a canonical momentum 

P = d^S. (18) 


Since p is not the classical canonical momentum, but an average, the theory 
remains fully quantum mechanical [24) . Particularly, the Hamiltonian 

dS rP 

T^Bohm =-= T. - \-V~\-Q (19) 

at 2m 

has a quantum potential Q Emn]: 


Q = - 


hP dl IV'I 


2m I'i/'l 

The form ([Tn|) can be re-expressed (using i? = |^| = as: 

/ da;R'' 

A dx \ p J 2 dx \ R 


( 20 ) 


Vp (j)\x) = -■ 


( 21 ) 
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The imaginary part of the weak value m is 

. r . 1 _ _ d^R 

Im [xP^] — 2 p ~ ^ R ^ 

where the experimental p term from @ is suppressed. Equation (ED then 
indicates that the weak variance is proportional (with a dimensionful factor 
of fi/2) to the divergence of this imaginary part. 

Applying the divergence operator in (I^Tl) gives 

and by recognising each term above (with (l20ll and (l22ll l. we find 

Vp(pk) = f^^^^+mQ. (24) 

This nontrivially relates the weak variance to the quantum potential of 

the Hamilton-Jacobi picture, providing further evidence of the physical im¬ 
portance of Vp. The imaginary part of the weak value also appears in the 
expression above. By combining (ED, (EID, and (ED we obtain a nonlinear 
differential (Riccati) equation fixing the form of luixP^ at any given time, 

dlm^p^ 1 2 2m 

=-j-OW (25) 

This will be an important element in our discussion of the classical limit 

(SID. 



2.1.2 An expression for Vp {p\x) in terms of weak values 


Consider the weak value of momentum squared. Substituting once again the 
polar form if = |'0|e®'^ and gives 

Recognising hdxS as the real part of the weak value of momentum, we can 
re-express the quantum potential Q in a well-known variance-like form (up 
to a scaling by 2m): 




{x\p'^ \tp) 

{x\iIj) 


Re 


{xIpIiP) 

{x\fj) 


( 27 ) 


Now consider the weak value of momentum, squared: 


Re 


(xlpl-fp) 

{x\i}) 


= Re{xPp)'^ - ImixPp)"^ = {hdxSY - Im(a;p^)^. 
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Clearly, we can subtract this from (IMl) above to obtain 

(28) 
(29) 

by use of (IMl) and linearity of Re( - • •). This expression depends on both 
the first order {xPip) and second order {xP^) weak values, and resolves the 
ordering ambiguity, discussed in the Introduction, between taking the real 
part and taking the square of the first order weak value (compare (1271) and 

(IMl)). 


2 mQ + 


lTCL{xP^f 


= Re 


Vp {p\x) = -Re 


{x\f H) 

{x\'ip) 
{x\p^ \i>) 

lx\ll)\ 


— Re 



2.1.3 The weak variance is not a weak value 

We have reached an expression for Vp {p\x) which superficially has the form 
(p^) — (p)^- Let it be noted, however, that this is not a construction in terms 
of moments of a distribution: the complex weak value xPip is not a moment, 
it is the real part p = He{xPtp) that plays this role. 

Consider, now, expanding the square in (ITT]) to obtain 

'^p(pk) = (^ Jp'^yV{p\x)dp^ - ^pf =Ep(/|a;)-E(p)^ 

This also almost looks like (p^) — (p)^, however the two expectations belong 
to different (pseudo)distributions (>V(p|a;) and yV{x,p) respectively). Using 
the law of total expectation, we could write both terms using Epi 

Vp {p\x) = Ep(p2|x) - Ep(Ep(y|p))^ 

where y\p is (formally) a position postselected on momentum. This averaging 
over all positions y shows that despite being locally valued (in terms of the 
postselection x), Vp {p\x) encompasses information from the entire configu¬ 
ration space. This is not surprising if one remembers that the wavefunction 
solves a Schrbdinger equation over this same space. 

Quantum mechanical expectations are calculated in the Wigner formalism 
as integrals over the entire phase space, 

(A) = J dx J dpA{x,p)W{x,p), (30) 

where A{x,p) is the classical observable and Weyl operator ordering is as¬ 
sumed M- One may then define configuration space densities 



dpA{x,p)yV{x,p) 


(31) 
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which give the expectation {A) by direct integration, and ‘local densities’ [2 


^loc - IV’I ^cs 


dpA{x,p)yV{p\x) = Re 


{x\^) j’ 


(32) 


which give the expectation by weighting the integration by the probabil¬ 
ity distribution |'0P- By writing our classical observable as a polynomial 
A{x,p) = Yhn fn{x)p'^^ we find that these local densities are nothing more 
than linear combinations of conditional moments / weak values of p". We 
have already seen (1291) that the weak variance is quadratic in the weak value 
a;Py,, therefore it is not a local density and no classical observable associated 
to it exists. Furthermore this suggests (since classical observables are quan¬ 
tum observables too) that there is no weak variance operator Vp, as argued 
on interpretational grounds in HI. 3. 21 Therefore, the weak variance is not a 
weak value. 


2.2 Hydrodynamical approach to the Weak Variance 


The real part p of the weak value is the (Euler-picture) momentum of the 
Madelung fluid at position x. Classically, the fluid momentum is the (local) 
average momentum of a particle in the fluid; quantum mechanically, p retains 
this qualitative feature of being a local average as discussed above. 

In the Hydrodynamical picture, one can derive wm a pressure tensor 
using the same logic m from which one derives the fluid momentum (HZl). 
We obtain in one dimension (cf. [2], Appendix C): 

P = — [ {P - pfyV{x,p)dp = — (-^d^ln{p)) . (33) 

m J m \ i J 


Combining the pressure above and the equation of state of the Madelung 
fluid P = pksT gives a definition for the thermal energy [2: 


keT 


k? 2 

— -—dl Inp = 
4m 


Vp {p\x) 


(34) 


There is an intuitive picture relating fluid temperatures to momentum 
variances: In the kinetic theory of gases, a higher temperature corresponds 
to a wider Maxwell-Boltzmann distribution of momenta, according to cr^ oc 
rnksT. The corresponding quantum relation (13411 preserves this intuition at 
the fluid level, for a temperature distribution T = T (x) over postselection x. 

This expression dSSD should be compared with the independently mo¬ 
tivated expression for the weak variance. In this light, the equation 
of state may be justified (rather than postulated) by writing the Wigner 
function from dSSl) in terms of its marginal, using m to find the solution 
P = ^Vp (p|a;), and then invoking this analogy to kinetic theory to link P, 
p and T by an equation of state. 

Physical potentials tend to arise from interactions (mediated by gauge 
fields in the Standard Model or the metric field in General Relativity), but 
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the quantum ‘potential’ Q puzzlingly appears to pop out of the Schrodinger 
equation. Interpreting the weak variance as a thermal quantity pushes the 
‘interactions’ which give rise to Q into a hypothetical microstate description 
from which postselected quantum mechanics emerges thermodynamically, 
potentially resolving this philosophical worry. 

We purposefully do not propose at this time any underlying statistical 
thermodynamical picture, instead pointing the reader to the discussion of 
’’subquantum theory” in [5], which suggests that one should analyse the 
properties of the Wigner Function (such as, in this study, its weak variance) 
in order to gather some information about such a theory (cf. >13.41) . It is a 
strength of thermodynamics that it does not formally require any model for 
microstates. 


2.3 Measurement of the Weak Variance 

We have shown that the weak variance can be described without reference 
to the properties of a specific apparatus, using only the description of the 
system’s state through the Wigner quasidistribution or the wavefunction ip = 
. However, combining the many equivalent expressions developed in this 
section with the t] parameterisation of the imaginary part ([2]), shows that 
this quantity should in fact be easy to measure in a number of experimental 
settings. 

Unsurprisingly, the weak variance of momentum can be measured using 
weak measurements of momentum. The form ([291) is available if one has 
access to the second-order weak value xP^- With a good calibration of the 
apparatus’ rj and a good resolution over x, (EH can also be used. 

More interestingly, weak measurements of momentum are not necessary 
to quantify the weak momentum variance. Given the probability density p (x) 
determined from a series of position measurements, one may simply use the 
—dx In (p) form to find the weak variance at any position, without performing 
a single momentum measurement, weak or otherwise. A combination of weak 
momentum measurements and projective position measurements (possibly 
the same measurements used for postselection) can be used to cross-check 
the previous two measurements using a single dataset. 

One might also verify experimentally the law of total variance (1151) : for 
example, one can measure the total variance V (p) from projective momentum 
measurements, and one can measure each of the partial variances from weak 
momentum measurements. 

Having these many experimentally independent ways to measure Vp {p\x) 
should allow this formalism to be stringently tested by contemporary instru¬ 
ments (modulo a calibration or model for the apparatus parameter p ). 
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3 Negative Weak Variances 

Quantum Mechanics predicts that in some regions, Vp {p\x) < 0; this odd 
behaviour is known to the literature, but often glossed over (cf. e.g. the 
footnote to the Appendices of our primary source m)- The Bohmian and 
Hydrodynamical formalisms even allow the ‘weak trajectories’ x (t) to cross 
from regions where Vp {p\x) is positive to regions where it is negative. It is 
not sufficient to say ‘weak values can lie outside of the eigenspectrum’ to 
explain why a variance can turn out to be negative. 

Furthermore, this is not just an odd but unmeasurable feature of the 
theory, like the negative probabilities of realisations of p|a: [T UTTlfT^ . The 
second cumulant may, as we have just argued, be measurably negative. 


3.1 Explaining Vp {p\x) < 0 

The reason we identify Vp {p\x) as a (conditional) variance is that it ap¬ 
pears as the second term in the expansion of the (conditional) characteristic 
function into cumulants (|8]). We therefore consider to what extent such an 
interpretation for this quantity is warranted beyond the (experimentally suc¬ 
cessful m) identification of the first cumulant as a conditional average. 

The conditional variance is given by the second central conditional mo¬ 
ment of W (ED: 


ip\x) = y ^ ~ ^p^ ( 35 ) 

which may be contrasted with the marginal moment (ED- The term {p\x — p)'^ 
is the square of a real quantity, so it must be positive. Any negativity in 
Vp {p\x) must therefore come from the distribution >V {p\x)] and indeed, the 
conditional Wigner function can take negative values [1]. This should be 
contrasted with Yx (p), which depends on the marginal distribution W(x) = 
> 0 and thus remains positive. 


3.2 Wigner Negativity and Standard Deviations 

Even though the conditional variance may be negative, allowing imaginary 
standard deviations of p\x would violate the requirement that p,x gM. (i.e. 
the requirement that realisations of a real random variable cannot be com¬ 
plex). A standard deviation defined using an ad hoc extension to negative 
values, 

ap\x = \Jahs {Yp{p\x)), (36) 

satisfies ap\x G such that all our random variables are real, as required. 
This extension manifestly reproduces the usual standard deviation when 
Vp {p\x) > 0. This last point is necessary not only for mathematical con¬ 
sistency, but also for consistency of the thermodynamical picture (13411 : The 
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standard deviations of momentum in kinetic theory are obtained from the 
Maxwell-Boltzmann distribution a oc y/mkBT. 

However, the conditional distribution is a pseudodistribution, invalidat¬ 
ing the simple expression (I5S1) . Consider the integral form (1551) of the weak 
variance: even when Vp (p|a:) = 0, there can be a spread of p\x away from the 
weak value, since the negative and positive parts of the pseudodistribution 
W {p\x) may cancel under the integral defining a cumulative distribution, 
without any information given about the squared-term (p\x — p)^. Expressed 
more formally, deviation-bounding theorems such as the Chebyschev inequal¬ 
ity do not apply to quasidistributions. This possibility for stochasticity at zero 
variance is further emphasised by the fact that higher n > 2 cumulants (EH) 
are generally nonzero, indicating that the pseudoprobability distribution is 
generally nontrivial. 

One might then attempt to push the absolute value inside the integral that 
defines Vp {p\x) to take care of this problem, defining instead our standard 
deviation with 


= J abs(>V ip\x)) {p\x - Ep {p\x)f dp. (37) 

By virtue of the triangle inequality |/ Adx\ < j \A\dx, this standard 
deviation gives at least as large a spread in p\x as Furthermore, it 
clearly reproduces (IMll whenever W {p\x) > 0, and so inherits all the desirable 
properties discussed above. 

However, this approach is similarly unacceptable. The expectation Ep is 
also an integral over p, yet if we want to consider deviations away from p then 
we cannot replace the conditional distribution by its absolute value in this 
expectation integral: changing the distribution changes the random variable 
of which we consider the mean! Furthermore, the correct normalisation of 
the quasidistribution J dpW (p\x) = 1 guarantees that abs(>V (p\x)) is not 
normalised, and so it is not a proper probability distribution either. 

In light of these failures, it is even possible to dispute whether a standard 
deviation is a sensible quantity for quasidistributions. Unlike the variance, 
which pertains to the distribution, the standard deviation is a deviation of 
realisations of a random variable from their mean. As discussed previously, 
realisations of p\x can occur with negative probability: building a sample 
standard deviation from these phantasmagoric realisations would be mean¬ 
ingless as the final step of a calculation, even in a quasistatistical outlook. 


3.3 Shape of '0 

In the description of Bose-Einstein superfluids, a quantum-potential-like term 
appears in the polar decomposition of the Gross-Pitayevski equation. In this 
context, Q is interpreted as the kinetic energy cost of curvature of the con¬ 
densate wavefunction - it costs energy to squeeze a wavepacket. 

Given that the weak variance depends in part on the quantum poten¬ 
tial (l24l) . one might wonder whether (by analogy) Vp (p\x) is related to the 
curvature of the wavefunction. 
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We now consider how the curvature of the wavefunction amplitude, R = 
I'i/'l, determines the sign of the weak variance. The relevant quantity is the 
convexity d^R of the amplitude, which is positive for convex intervals of |'(/'| 
and negative for concave intervals. 

3.3.1 Non-Nodal Points 

We use the decomposition (l23ll to solve the inequality Vp {p\x) > 0. We note 
that R is nonnegative by definition: for non-nodal points (i.e. points x such 
that R{x) ^ 0) the weak variance is nonnegative whenever the convexity 
satisfies 

> dlR. (38) 

The left-hand side of (IMl) is positive, so concave regions (0 > d^R) always 
have positive weak variances (independently of the system under considera¬ 
tion) . Regions that have a very small concavity are also very likely to satisfy 

dMl- 

Maxima of the wavefunction amplitude are concave so Yp {p\x) > 0. Near 
local minima of R{x), the left-hand side of (IH^ vanishes and the inequality 
is never satisfied (since d^R > 0): Hence we expect Yp {p\x) < 0 near local 
minima (again, independently of the specific system under consideration). 

In the extremal case, piecewise exponential amplitudes R (x) = exp(±fca;) 
always give Yp {p\x) = 0. These amplitudes are observed e.g. for particles 
tunnelling into a constant potential barrier. 

Since the position-dependence of the amplitude R{x) is easy to measure 
with projective position measurements alone, this inequality can serve as a 
test for the sign of Vp {p\x) even in an experiment not optimised to calculate 
the weak variance. We recommend that experimental probes into negative 
variances should occur near non-nodal minima of the probability density, 
since the convexity in other regions may not support negative weak variances. 

3.3.2 Nodal Points 

A slight complication arises around nodal points R —>■ 0 of the wavefunc¬ 
tion: The weak variance at the node itself is infinite, given the logarithmic 
singularity at ln(0) that appears in (ITU)) . 

To show this, we note that continuous differentiability of ip and ijj* en¬ 
forces dxR —>■ 0 around nodes. The expansion of the amplitude R about a 
node at any fixed time is then 

R{x) = ax'^ -I- 0{x^), 

where the condition > 0 (that a node is a minimum of the amplitude) 
gives a > 0. Around this node, the weak variance (ITUl) tends to: 

— ^d'^lnax'^ = —h^d^lnx = —>• -|-oo. (39) 

It is simple to see that this asymptotic result generalises even to pathological 
expansions i? = 0-I-- • --l-O-f 0(a;"'). Interestingly, this behaviour for the weak 
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variance is always strictly positive, excluding Wp {p\x) < 0 very close to the 
node. This will be illustrated in the following section in a number of systems. 


3.3.3 Illustration in toy systems 


The study of coherent states in the context of negativity of the Wigner func¬ 
tion is motivated by the result |28] that only these states are nonnegative 
everywhere. 

Consider the coherent states of the Quantum Harmonic Oscillator. We 
must check (for consistency) that Vp {p\x) > 0 (strictly) in these states. The 
weak variance is 


Vp (pIx) 



and, after many cancellations, one obtains the simple expression 


hu! 

Vp [p\x) = +m—. 


This is indeed positive for all x, as expected, and manifestly has the 
required dimensions for a momentum variance: 

[mass] X [energy] = [momentum]^ . 


We note that Vp {p\x) jm is equal to the zero-point energy Sw/2 in this 
minimum-uncertainty state. In fact, for any energy eigenstate of the harmonic 
oscillator, we hnd a zero-point contribution to the variance: 


Vp(p|a^) 

m 


hu) 

'Y 



InHn {y) 



(40) 


and this contribution remains even for arbitrary superpositions thereof. 

For an energy eigenfunction of a particle in a box (respectively in a — a/ |x| 
potential), we obtain 


Vp {p\x) 
m 

(resp.) 


+Er, 

2m 


2 , T17T , -, 

CSC (-^F - ^oj) 


■ Er, 


IuLI^^Ii iy) 


(41) 

(42) 


in terms of the energy levels En,n > 1 of these two systems. The eigenstate 
in a box or a harmonic oscillator is clearly not in a coherent state; yet we 
find that Vp (p|x) > 0 for all x for these toy models, because each x lies in 
the vicinity of a maximum or a nodal point. 

We note in passing that the weak variance becomes positively infinite 
at the nodal points of these systems, as expected from the discussion in 
the previous section. More interestingly, the weak variance tends (at certain 
postselections x) to the zero-point energy when the system is in its ground 
state. We purposefully do not propose at this time an interpretation of this 
finding. 
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3.4 Negative Temperatures 

If the weak variance is the temperature distribution (IMl) of the Madelung 
fluid, then the Madelung fluid can have negative temperatures. Negative 
temperatures may be surprising, although their study is theoretically well- 
established [53] and the concept has successfully been applied to several ex¬ 
perimental systems exhibiting population inversion. Some interesting prereq¬ 
uisites for negative temperatures are that the internal energy of the system 
must be bounded from above as well as below, and that T = 0 requires the 
energy to be extremised to one of these bounds. The ideal gas law P = pksT 
shows that the pressure also becomes negative, precisely when the tempera¬ 
ture is negative. 

One further consequence of weak variances on the thermodynamics of the 
hydrodynamical formalism is interpretational. Recall that the hydrodynamic 
momentum p = E(p|a;) is interpreted [5] as the local average momentum of 
the fluid particles at a position x. Yet, we have seen 1 93.21) that negative weak 
variances do not provide a standard deviation or a Chebyschev inequality 
for the random variable p\x. In the hydrodynamic language, the peculiar 
momentum of a Madelung fluid particle is not required to be zero, even at 
zero temperature! Also note that the temperature of a nodal point in the 
fluid (where we expect to find no fluid particles) is -|-oo. 

These are the sort of features any “subquantum theory” [5] should exhibit, 
if such a theory admits the Madelung hydrodynamics as its thermodynamic 
limit. 


4 Nondeterminism in Quantum Mechanics 

We have already discussed ( 93.2|) how a null weak variance still allows stochas- 
ticity and nonzero standard deviations. Although Vp {p\x) = 0 does not imply 
determinism, one might nonetheless consider whether a relation between the 
weak variance and determinism exists. For example, the highly nonlinear 
form of the weak variance makes it unlikely that a null weak variance would 
emerge for superposition states (or, a fortiori, for entangled states). These 
conditions correspond exactly to the nonclassical regime, where we expect 
nondeterminism. 


4.1 The Classical Limit 

4.1.1 Consistency of the Classical Limit 

It should not be surprising to find that these worrisome negative variances 
disappear in the classical limit. Heuristically considering the classical limit as 
the limit of infinitely thin coherent states similarly guarantees Vp {p\x) > 0 
since the conditional distribution is nonnegative |28] . 

However, a stronger statement is also true: weak variances should dis¬ 
appear completely in the classical limit, since classical observables are well- 
defined and motion is deterministic. We have even shown f 92.1.3|) that there 
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is no classical observable associated to it. The quantum potential behaves 
like Q —?> 0 in the classical limit 0- When <3 = 0, the Riccati equation (1^^ 
admits a trivial solution IvtixPip = 0 since the weak momentum (g C) tends 
to the classical momentum (g M). Hence, using ( 1 ^ . Vp (p\x) — >■ 0 classically. 


4-. 1.2 Deriving the Classical Limit 


Another approach to the classical limit would be to minimise the contribution 
of the scedastic function Vp {p\x) in the total variance ((T^ by our choice of 
wavefunction, in order to mimic determinism. Since the weak variance does 
not depend on the wavefunction phase, only the probability density p = 
is relevant. We can express this extremisation problem in the language of 
calculus of variations as 


6J[p]=SiE, iYp{p\x))) 

= 6 f dxp{x) dl \og{p{x)) 

(.y 


= S dx 


P 


= J dx S£{p,p',p") 


= 0 , 


(43) 


where unimportant constant factors are suppressed. The resulting Euler- 
Lagrange equation depends (albeit trivially given our specific C) on the sec¬ 
ond derivative p" = d'^p of the probability density: 


d£ p" 1 (p'V 

dp dx dp' dx"^ dp" p 2\p) 


(44) 


We notice that the terms of (1231) are similar to those of the Euler-Lagrange 
equation (PI). Consequently we find that in this limit, a term of our choice 
cancels out of (1^ to give 


Vp {p\x) = (Im [xPp\f = -mQ. 

However, reinserting this into (P)) results in Vp (p\x) = -\-mQ] this is uniquely 
possible for Vp {p\x) = Q = Im [xPp] = 0. We have therefore derived that the 
minimum average contribution of the weak variance to the total variance is 
zero: negative weak variances do not allow the nonweak variance 'Vx{p) to 
exceed the total variance in m- 

In summary, by minimising the stochastic part of the total variance we 
can derive the classical limit Q —>■ 0 of the Hamilton-Jacobi formalism. 

^ One must strictly take the limit of the quantu m p otential for the entangled 
state of the system and the measurement apparatus |22| , but Q ^ 0 for the system 
alone is not seen as problematic in the literature. 
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4.2 Vp {p\x) —)> 0 as a Semiclassical Limit 

Equation (EU states that Vp {p\x) is the divergence of the imaginary part 
of the weak value: when this variance is zero, Im( 2 ,p^) is a function of 
time only. In this case, either (I24p or the Ricatti equation (051) . gives that 
Q{x,t) = (5(tonly) < 0. Hence the quantum force dQjdx = 0: this is [53] 
the condition that Bohmian trajectories coincide with their classical counter¬ 
parts. Vp (p|a;) —>■ 0 therefore yields classical particle dynamics (on average), 
occurring at nonclassical energies. Of course, these nonclassical dynamics are 
localised to specific postselections a;: if Vp {p\x) = 0, Vx then the average 
weak variance is also zero and we recover the classical limit as derived above. 

Consider a particle tunnelling through a constant barrier: The exponen¬ 
tial wavefunction amplitude guarantees that Vp {p\x) = 0 inside the barrier, 
therefore the tunnelling particle lives in this semiclassical limit and will (on 
average) follow the path of a classical particle with a higher energy: it will 
go over the barrier. The fact that tunnelling particles behave semiclassically 
is indeed why tunnelling problems are successfully visualised and formulated 
in terms of particles. 


4.3 Weak variances do not ‘measure’ a deterministic particle ontology 

It is sometimes claimed that the weak value Re(x_Pi/i) ‘measures’ the Bohmian 
momentum [Bill]- A naive Bohmian might expect that the weak variance 
would vanish, being the variance of the nonstochastic Bohmian momenta. In 
this section we argue against one possible misinterpretation of these quanti¬ 
ties: weak values (and their associated weak variances) do not measure the 
properties of a Bohmian particle. 

The only hidden variable of a Bohmian particle is position - it does 
not even have an energy or even a momentum independently of what the 
mathematical formalism predicts we should measure |221l30j . This is of course 
a reflection of the fact that nonposition measurements in Bohmian mechanics 
are contextual [22] . As such, a measurement of the weak value, although equal 
to the Bohmian momentum dxS, is not a measurement of the “particle’s 
momentum”; similarly, the weak variance does not measure the standard 
deviation of Bohmian particles’ momentum (both on interpretational and 
quasistatistical grounds). 

With the above caveat in mind, the weak variance allows a determinis¬ 
tic particle interpretation of the ‘semiclassical’ limit Vp {p\x) —>■ 0: such a 
(localised) Bohmian interpretation has already been shown to apply at post¬ 
selection points X where it is at its strongest (particles tunnelling and in the 
classical limit), and is guaranteed to not apply where the global Bohmian in¬ 
terpretation is known to break down | 551 [ 531 [ 3 n[ 35 ] (since Vp {p\x) —>■ -boo ^ 0 
at nodal points). We emphasise that a semiclassically emergent Bohmian on¬ 
tology is neither supported nor precluded by the absence of a sensible notion 
of standard deviation for quasidistributions, nor can it be of any instrumental 
relevance. 
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Conclusions 

The weak variance can be defined naturally and uniquely in terms of the 
Wigner function, singled out by insisting on a (quasi)statistical interpreta¬ 
tion of quantum phase space. A weak variance does not define the standard 
deviation of individual realisations of p\x, which may occur with negative 
probability, but instead describes the system’s state. 

The weak variance Yp (p|a;) contributes to the momentum variance of 
the uncertainty principle, justifies the ‘ideal gas’ equation of state of the 
hydrodynamical formalism, and appears to be related to the ground state 
energy of a few analytically soluble systems. An expression Ke^xP^ — xPip'^) 
in terms of weak values was derived, which we used to argue that the weak 
variance is not itself a weak value. 

We show that the sign of the weak variance is predetermined for local 
extrema of the probability density: Yp {p\x) is positive for nodal points and 
maxima, and negative for non-nodal extrema. We also derive some features 
that any hypothetical ” subquantum theory” (from which quantum mechanics 
emerges as a thermodynamic limit) must satisfy. 

We show that the weak variance is not only consistent with the classical 
limit, but can be used to derive the classical limit from a variational principle 
where the expectation of the weak variance controls the emergence of deter¬ 
minism. This construction also shows that the variability Yx {p) of the weak 
value cannot exceed the total variance of the state. A semiclassical limit is 
derived and discussed. 

Most importantly, we note that Yp {p\x) should be measurable by a num¬ 
ber of experimental setups (' ^2.31) . 
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Appendix: Derivation 


The characteristic function of a distribution is its Fourier transform. For the 
Wigner function we can write the conditional characteristic function [mn] 
as follows: 


M[t\x) 


J dpyV{p\x)e^'^P/^ = 


^*ix - + ^t) 

'ip*4’ix) 


(45) 


The cumulant-generating function is the logarithm of the characteristic func¬ 
tion. This gives 

ln(M(T|a;)) = ln( '^ + ln('0(a:''')) — ln(p(x)) (46) 

p[x) 

= ^ [Mp(a^“)) + ln(p(a;+))] -b i [5'(x+) - ^(x")] - ln(p(x)), 

(47) 

which we can then Taylor expand to obtain the cumulants [T]. 

By inspecting m. we expect final factors of {hj2i)'^ and an expansion 
in /nl, so we consider series of the form: 


h 


OO 


/(a^± ^(*t)) = ^ 

n —0 *- 


±- 


'2i 


(d"/)(x) 


(iry 


where the (fi/2*)” come from the chain rule. 

From the i{S~^ — S~) term we obtain a contribution to cumulants of 


i 




^d^{hS) if n odd 
0 if n even 


and from the In p terms we obtain: 


1 

2 




d” Inp 


'O ifn = 0 

< (H) ” d" In p if n even 
0 if n odd 













